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Description 

The p^^nijnv^nton relate so computer^ed tomography scanning and in p&t$culm t to an improved method and 
apparatus for&cajniraog s subject and constructing a visual image from the .sensed tomography data. 

Garrpui&iz&j lomography has been used for many years. particularly m tie medical fte&, providing * 
s&uotsve or ftm\r&&si*(s method of gensmtiog visual sectional views of a Subject such as a patent, as described m U,$ 
Patent No 5.414,62$?. the entire contents of which is Saccsrporatsd herein by reference. These object scanning devices 
process semed tomographic data and provide a visual image of the subject as illustrated by US, Patsm Mo&. 
5,398,684: 5,099,346; 5,391,877; 5,414,628 and 5,402,462. 

It is an object of the peseni invention to provide an improved conpuier^ed tomography scanner and an improved 
method of coo^putehzsd iooxK^sphy scanning. 

Another object of the present invention- is. to provide an improved computerised tomography scanner and a method 
of conputehsed tomography scanning to produce Images of better quality than conventional scanned. 

Mother object cf the preset Invention is to provide an improved method and apparatus incorporating a transform 
inversion scheme and a reconstruction slgoothm of redded computational complexly 

H is further an object of the present invention to iitiU^e a multlscale backprojt^lon with postprocessing step instead 
of a conventional backproleo^on algorithm or direct Fourier method to obtain Improved images. 

Vet another object of the present invention is to provide multiset methods wh^ can be appiied under weaker reg- 
ularity requiremer^ than Fourier methods, so the present algorithm can be adjusted so provide different resolution for 
different parts ol the reconstructiof^ whether or not the Radon data are equally spaced. 

With these and other objects and advantages o? the pi*esenf invention as may become herein apparent, the inven- 
tion may he more dearly understood ffom the following description, the attached figures, and the upended daims 
herein. 

m\BF DESCRIPTION OF THE DB.AWlMClS 

Fig 1 illustrates the Radon transform cf f at angle- 0 and position % denoted IVit), which is the Integral o! i on Ihe 
line x cos § * y sin a « t . 

Fig, 1 A sllusirates is a hardware block diagram in accordance with the present Invention. 

Fjg. 2 shows a construction of sero-level grks g 9 « whloh contains the results of a single bacteroieclion. that cf 
according £o a prof erred emSxxSmeni of the present Inventor 

Fig 3 shows the merge of and g.^ k to form n acooming to a preferred embodiment of the present inven- 
tion. 

Fig, 4 shows a choice of spacing in the slow direction according to a piBferred embodimem of the p? esen* invention. 

Fig, o shows the &$ftH8& symmetry of the Radon transform - convolution bac&projection sequence according to 
a preferred embodiment of fhe present invention, 

Fig, 8 shows op&nal values a* the Gaussian width according to a referred embodiment of the pnssent invention. 

Fig. 7 shows surface plots of point spread functions for the classical and mulfiscaie bacl<projection methods and 
surface plots cf their Fourier transforms. 

Fig. a shows an Image from which Radon data were computed for the purpose of testing different variations of com 
volution baeknrojectlon (CB). 

Fig. 9 shows the results of CB with classics backpro|ection. 

Fig. 10 shows the results of 08 w3h mu&scale Pad^projecflon without postprocessing. 

Fig. 1 1 shows the results of OB with mu&scaie bac^rojeotson, without postprocessing but with doling of the data 
by means ol interpolation. 

Fig. 12 shows the results of OS with mulilscaie baci^ ejection, without doubting of the data but with postprocessing 
accenting to a preferred embodiment of the present invention. 

Over the past decade, numerous atieoipie have beers made to improve the quality and speed of various techniques 
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lor oon^uterfe^d lombgr^y scenting or othsr object image scanning Such t^ohniques are no* mmed strictly to the 
mdrffcfei ftek*, not may fir*d application in other areas such as synthetic aperture nad&r image reconstruction in meny 
syndetic aperture radar systems, ssmipiea of the 2-0 Fourier transform of a radar r^iecMv% fynctloo are obtamed on 
a polar grid in tha Fouher plane, $nd backprojecfson may he used to reconstruct the reflectivity function from these $am* 
« pies. 

Condenses fomogfapfty scanning Is a nondest? active or noninvasive method at generation actional of ao 
object. Computerised' karKxp^y or CT scanning has bean appitad In induatnal processes tot flaw detection and other 
analyses -and baa bean adopted by the medical community *o? the examination of patients. 

Figure 1 A is a diagrammatic representation or the computerized tomography scanner of she present invention. Tha 

re object s&anner 2 includes a stationary frame lor mounting a source of eleof ronvacmetlc radiation ■ sach as- X-rays, to scan 
a subject and an array of detectors or sensors positioned to receive the scanning radiation. The radiation source may 
produce a l&n-shaped beam which la directed through' .an object to be examined towards the detectors. The defectors 
of recede radiation fen the scarce and produce data in tha form &f electric signals indicative of me intensity 

of the mdiafion Because the intensity of the mdlafson from the source is known, the intensity of tha radiation Shrptnging 

ts ypon en tntiMtiusi sensor is representative of the attenuation of the radiation beam by the ablest along the path 
between the radiation source and the detector T he eseetheai signals indicative ei this Intensity from each of the Individ- 
ual detected or censors are subsequently conveyed k> a promotion data memory 4 for storage. 

The scored signals are inputted to an Image processor 6 ior processing the representation of tha raolatton attend- 
atien aom the Individual sensors. A caicaiating device 14 within the insane processor computes a mathx of samples 

$o from the s^nais supplied by tha sensors or detectors. The processing means provides a matrix of samples of the image 
which are stored in the image nsatnx storage 8. Each detector sampling represents attenuated data collected aiong one 
of the paths of the projection from the radiation source The image matrix storage element & connected k> a video mon- 
itor IS and pdnsnp; device 10 so that a visual Image can be produced. 

The Radon Transform m a known mathematical operation which is generally equivalents computing the integral of 

s& a function along a given Sine and has bean widely used In computerised tomography Thus, the inversion of tha opera- 
tion Is parttetiiafiy useful In making radiographs of plane sections of a body or object. The pus pose of tomography is to 
show detail in a predelera^ned plane of the body. Reconstruction of & function of two or three vsnabses from *rs Radon 
transom has proven vm in computed tomography (CT), riutiea? magnetic resonance imaging, astronomy, geophys- 
ics, and a number of other fields. S,R Deans, The Radon n^m^rm and Some of m Appfiomfom.Hmt York: John 

as Way and Sons, 1 983 One of the best known reconstruction algorithms is the convolution baol^mieciieo method (08).. 
which :S widely used in commercial CT devices (with running for divergent^esm projections): B. Hmtmn, im^ge 
BecornrtriKibn from Prohct&ns: The-Fhndammtakt of Computerised Ibsn^n^phy, Mew York: Acade?nk; Press,, 1980). 
Recantly : it has been appfced to spot%hFmods synthetic aperture radar image reconstruction, In which the conven- 
tional method Is the direct Fourier method (OF). u&.., Fourier domain interpolation flowed by 2 -D FFT 

@g Convolution Backpro^ction (CB) Is preferred to DF since the termer provides better images. R- Mersereau end A. 
Oppenhelm, "Digital racons?roc^oc of n^ltldirnenslonal signals from iheir pjojections," Pmcea^np? o^ the lESB, voh 
62, no, 10. Oct.. 1974, pp. 1319-13^8. However, smce the bad<prdi^t?on part of CB raises the con^atatessi complexity 
of the method to 0(f#}.. wtifie DPs complexity Is 0(M s log M), there has been interest in finding effective implementation 
o^DF 

By way of fdrther backpmund, there ere a nunt>er of mathematbal considerations which should be conslde? ed with 
re^>ecl to the invention described below 

For the purpose of discussion, let f be an absolutely Imegrable function of two verges, P$j$< the R^on $rans- 
fern ot the fuhofon at angle a. Is then defined by 

F 0 f( i } f f{t cos a - x sin a, i sin a e r cos S.) m-, CO 



m This is the Mite^ra^ of f along the iine 

x'cos 0 + y sin 0 » t 

in the x -y plane. (See Fig. n, \n many Imaging technoio^es, such as CT, the challenge & to compute f, or an approx^ 
matlon to it. given sample of P^^t). for finitely many values of n and t. 

To understand c:B. which ;-s meant to solve this pmhiem, one needs the Fourier tmndpnn. For f v. C{ E °). , the Fou- 
rier transform f of Ms defined oy 



ESWSsXX^iD. -:CP 



A 1 

( 2-jt) w * ' «" 



for ®H u t: R d . For any a, / $ : I U M 2 ) imp»es that P$f is w&8 defffi&tf and belongs to L 1 (R v ), and thus 



rise Fcsuriar transbrm of ?V&) wlh f«$|}$ct to i tea&o weft defied. Many methods tof reeoj^rucft^a function from its 
Radon -transfers Ine&idsnq C& rafy or? the Radon Slice Theorem, which states that 



p) « v/2«f(pcos paid #) 14) 

Thus, it fC&y} «s the knage function, the 2-D Fourier transfer?)? of the smaga can be @3&pfed by sampling the Fourier 
trah$#3?m-of the Radon transtor m of tha image to* cMemnt vaiues oi p and 

One way to reconstruct f from its Radon transferer? is to compute r i$s*ng {£} for as many v§Uu@& of p and 0 as ps*> 
sirie and th^h p&dorm a 2~D invsrse Fourier transform. AJgonthros based on th& idea are sailed direct Founar methods. 
(OF). SB, however, uses tne1o$iow*ng approach. 1 he Fs^t^ transform toarston forrf&fa for functions nr 2 y&rk^ksscan 
be written 'in polar ooordshataa as 



By th«? Radon Slice Tneoram (2). tNs ea&afe 



1 f *T~aa f / n N i & 1 ft ***** * J ^ rnufo d® 



In. most appliss&orm vainas of P< ; r' or 




are obtained from th$ physical data. Uesng trso 



P f 



valuea : one may compute a ^sx^ion 



> f 
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and then 



which means that f (x : y} can be reconstructed the F $ j/ vaiues, giver? these valises for a& 

Of course, In. actual applications. P* 5 f{§ is kncwn for only a finite number of values of 0 andi and f is to be computed 
at a finite sat of points, or pixels, (xy). Los us assume thattfta pixels {x it ys). to be computed &e inside a circle inscribed 
in anNxN square gdd, with a distance d between adjacent pixel centers in both the horizontal and vertical directions, 

20 that there are Q evenly-spaced attgiss 0, at which Fy" is known, and that at each such &ngle ©$> P$ff$ Is known at N 
evenly spaced values of t tft& drrlersnee between consecutive values % and V, ; hein§ d as well, Limiting the region 
of destruction to the Inscribed circle, N samples of the Radon transform In any direction are enough to cover the 
region. This exposition assumes that O Is a power of 2, though the new algorithm can de&ny be g&rttrai&ed to the case 
or arbitrary Q. 

&s Computation of 

j 



is catted the convolution step because multiplying two functions together In the Fouler domain and- taking the hw&tm 
Fourier transkxm ofthe result is associated with convolution. The R>urier transform end mverse Fourier transom 
r$qus*feti in (3) for this step are approximated in practice using the discrete Founer transform (DFt) and its inverse. 
Since there are M values of PrJit) for each value of 0, the computation of 



p { 

4* for a single value of a Is D(N irxr N) using the Fast Fourier Transform (FFT), and since there are Q values of a, the entire 
convolution step is Q(Gh&&3N). ii as is -usually assumed, G* N, ^en the e^>vt>Mioh step is O<N*ie0hi}, 

The Integral with reject to 0 in (4) is reps&ced in practice by summation over all -available values of 0. Since there 
are OChP) olseis at wNbh s&rmnaikso is required: and Q values of ® : this final step o? oompming the right side of (4) is 
clearly G(GN*7, or, if Q * H 0<N S ). Since the contribution of P,/ to a pixsi <x,y) depends only cn the quantity 
x cos t**y sin 0 , P$H*$ ^ said to be :< smeared" or l^actooiected" > aionq the line x cos 0 + y sin 0 « u $$r each value 
of u. That Is, it is added to every olxel lying on that line Therefore, this final step is called the baokprojectlon step. 
Although the concept of bacferojectlon suggests an algorithm with a different loop structure 'than the one suggested by 
{4} - namely, (4) suggests an outer loop to go through ail the pixels and an inner loop to go through the values of 0, 
whereas the backpro|eot^on concepts subsets the revere - the hackprc*ed;lon concept does not change the. compu- 
te tafiona! complex Although DP is 0(.N^N). commercial CT scanners have traditionary used G8, despite its higher 
confutation^ complexity, in part because of the s^enor quality of the Images it produces. 

The present invention represents an improvement over conventional coniputeosad tan-spgraahy. scanners and 
methods for tomography scanning to provide an improved visual image with lower computational oornpiexity In the 
present apparatus and method, the Radon projection samples from the detectors are processed by the flowing steps: 



(a) filtering the samples of each Radcn pf^ecttoo sample vector r, where; « I^.^Q, t>y: 

I. .confuting the discrete Fourier transform r ! of r :i using an algorithm of order N iog M, wher» N Is the length of 
me vector; 
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ft munipiyi^g tbfc fih ataYHini of ? by J, lor j < W8. and multiplying the jth element of r by H-j, fas- 1 > I4& t where 
the elements of r s are numbered 0 through KM : and 

m. oon^tsogihe inverse discrete Fourier transform g <: ; of the. moc$&«d r* using an algorithm ofocdar N fog h?; 
■and thereafter 

<*b) pertxmln^ the klh level of for H « l. -Joc^Q, by compulse^ the ohd g) > for I *? 1,..,X&£ k , by 

yrKS q«*>} by means of rotation of coordinates, truer poreljon ; and addition as described In dotal? hereinefter 

The present Invention comprises a fast, mu&scaie baches ectlon method and rx>stprocesslng method for exaco - 
tion by an. electronic calcels&ng device Referring now in esters in Ihe fast muffcscate backprojaction method or trans • 
for^^convosi^on-bao^rc^Krilon method of the present iroeft&yr. 

The mu^evel backpr^estioh method relies on the following reasoning. In the standard b&oi<pre^ection algorithm, a 
single NxM grxi 5 e osed to sum the appropriate values of 

P f 

for j « 1,... A for each pixel. One can, ix^wever. start with Q ghds -sj? J * 1 ; . .-A and tor eseh j project F^f only onto 
g $ \ later, the Q brfferent grkte can be added together p*xe!w?ea to produce me final *mage. prosed can be used 
to" save computation in the lowing way. 

In what follow for a£ i and /, if w$ he a function of two cdntimious. variables x and .y, snb g* f will he the sat of 
campus of f j at a finite sat of discrate points m the x »y plana. With every grid r*j is associated an single &j . and the 
po^sln the V »y plane at which samples of f- are computed for storage m oj are located *rra> the lntemeot:ono of two 
sets of paralfe; fines, one set parage! to the vector (eosoj , sin Oj } and the Other set perpe?K#cu&r to It The grids g/\ 
fbr various /. will be eaiied the seroih-ievel grids, the grids g. ? f for various /. be called the firsHevel ^ide eto. 

Given two functions fi and /jj , and iftm respective sampie grids g$ and g/ , if fj' y * m mih& mm fi and 

, then the method of the invention compotes g m f * 1 , the grid or samples of fj* ? ; ov smerpoiatfon from- and g; and ffj, . 
Gonputation of g^!* ' celled a v grid merge". 

Po? /« i,. V Q ? set 

/?(jr, y) * /J /{Jr cos 6). 



I.e.. f c ; is ^efuFxSbn of two vasiabfes msuKing from the backproi^ion of 

P f 

eiong the Hnea x cos h * y eln 0 « Uor deferent values oft In the x-y plane. Cfeer^y fH does not vary along soon Jinee, 
erxl therefore, when co^iectin^ sampler of to form the grid ^7^ it & suf ident io ooj^ip^te one sanipfe and atcw It in 
the grsd tos- oaoh of the M values of t at which P^f(t} has bean cornputed, raiher than computing and storing N^- sample 
See Figure 2. As ohown in Fig. 2, although g^f is meant to contain the samples of the function f'^j. defined on a d&fcsn 
the piano, ar^d the disk contains OCN 2 } p^<eb"i sho^n in the figure as the intersections of def^d Ijnes, it is eo^io^ens to 
compote and store In g ( i sanp^as of f H taken at oniy 0{N) polnts> showr^ in the figure as l^ack dots, s^nce f c j only 
varies In the direction {cos % In fact, Is merely a UD array containing all the competed values of 

and is oorrplled by the convolvers part of the CB algorithm. Thus, ^though the tea&^f&Ston begins with Q orids 
instead of one. the iolal nurnber .of samples stored Is O(MQ) f The distance between sample pokm of for eny j is d. 
ft^e distance hetv^eerj pt£t\t& where 
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if: SL£ScY*0$ed 

Eventual tsewous *n&al functions f° t na sat ha added together to formes fxv&\ ima^o. As ac&ition Is cornnutfa- 
|iv«f and associative, they may b$ added together m any ord^ and the wfhod of the iny^otson chooser an order that 

ic ;aduoea the nurnber of necees&ry computations. S:oce Q is svan by assumption the f'mi isvel of grid marges vvsil con- 
sist of adding togsith^r pslrs of functions / '? f f t such thai 0: and am close to each oth^n and thus, the direction in 
which f§ is constant ie dose to th^ direction jrs which is constant The sum of these two functions, A ; mil therefore 
vary slowly In a certain direction, and thus it will only be necessary to stoso aarrspies of f) from a nandfoi of widely- 
spaced points along eaeb line parallel to the "slow* direction. 

?*? Specifically it may be assumed without toss of geniality thai 0 ■< *>: < it for all j end that the 0; are ordered in such 
& way that 0: < &^ for at; j * 1 „ ,.,G - 1 . lor ail j « 1 ^ , >Q } M. 




Fes- all k from 1 to Qffi. one adds todoiher the functions ? and / 4 to obtain fl . Since for any ■ I, /'f does, not vary 
in the direction- (- sin cos 0 <)> it follows that f am f >?* vary slowly In the direction 



whore o I W< the average of O^^. j and Therefore, f % Marios si&tiy In the 



COS 



0 1) 



ns dif action, and when fbnt&*§ g | , It is: only necessary to storo samplas of f jt from a widely -spaced points ak>ng each 
line in that direction. The spacing between sample points in the K fasf direction, i.e,, 



I 
k 



sm 



01 h 
& ■ 



wi aga*n da d. The samples are computed from p ? f.s and fay Interpolation, 

A^hou^frthe flrsMevoi grids g \ each contain rrsorsi samples man any of the serothlevet grids g % Share are half 
as many f ftsMsvei. grids as xenoth^evel grids. afKi it will be sho^n that the totai number of sanies in any level of gads 

Bv similar reasoning, for £ 1 v Q/4, add / ^ I s end / 4 ? together to ctnain f 2 : ih such a wey that / 1 v^rlee slowly 
In the direction 



A/5 



c 



V 



v^here ml Is Ihe avamge of aod m> s , and con&equemly. when computing and stonng sanies of? ? to form t 
only necessary to etc^e samples too? a tew widely-spaced points in the d^eot?on 



JS 



x 



fl 



0 
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The spacing between sample points in the fast direction, l a. 

o > 

(cos #7, sin 075 

will again be ci f : | wlls not vary as slowly In Its "slow" directors as f varies in m %km" dirae&m. and thensfor« $'j 
ws» require rmr& sample points m Itsabw direction than gd/ requires in Ns slaw direction, ^evertbela&s. there are ftatf 
as many second^evel grids g ; j as there are fk&iAetf&l grids g 1 , and as w\\ be shown, the tofai number of point values 
that must be computed and stored for ml the grids at any one level of merges is O^N :s: VOM}- Samples of t H to be stored 
mg'< arecomputad fromg^f :< and gf ; by interpolation. 

Continuing .In this way, on© rn&y construct a sequence of leveis of functions and correlating grids. The tactions 
at the Hh Is^sl are constructed from pairs of functions at the 6*1 }--th tevsl In such a way thai each I -th level function var- 
ies slowly In a certain direction, ar^d thus the corresponding grid requires only a few sample points ?n that direction. At 
the {kx&vQHh i«l there b only one grid, and this grid represents the sum of &8 the original grips , fhai ;s. the sum 
of aB the hadsproleetlcns of the P&f. This grid Is therefore the resuming reconstruction, Since 0(N S ) operations are 
needed to build the grids at each level, the overall cost of the algorithm Is O(N^ogQ). 

This is shown In Figure -3 which Illustrates a merge of g 2 l< and to form The slow directions el f 2 ^ , 

and P^* are 

( sin ^/.. t , cos ih.i~V* (~ sm ^21^ cos 0 2k) * i%M ( ' m l*> cos 0$ I** 

respective;^ where Is the average % ? J, ; and §.>V Sample points are reps esented In tie diagrams by black dots.. 
g has more sample points- in Its slow direction than % 2 £.* and $ z f k do m theirs. The distance between sample points 
In the last direction, indicated in the diagrams by the distance belweeo adjacent solid lines, is always d. As explained 
above, the low eonputallooal coiripfexity of the algorithm departs 00 the judicious choice of sample point' spacing and 
hence, the m'te of sample points., in a grid s slow direction. This spacing is. chosen based on the toiiowlng reasoning. 

Since at every *tmg& level r two consecutive functions f ^1 s.r^ ?|l from the previous level are added to fern f 
I ... It Allows that 3? consecutive original functions t 9 were added, step by s*ep, to obtain f i . These 2 5 original functions 
are referred to as ti^e merged original f orations. Thi slow dilutions of any two conseoutsve original 'kindkra f % ar>d f 
{ % differ by vJQ radians, and It follows from this that the slow directed of Ihe f irsi and kst of the ?^erged orlginel func - 
tions differ by (2 : • l)?v/Q r&di&os, It is trivial to show thel the slow dir estion of ? c | Is halfway between the slow directions 
of the first and fast reread original functions, and thus ii differs by no more tf>an {2 { -1 } ;ri(2Q} from the slow direction 
of any of the merged ongina? functions. 

For any f i ^e last dir^bo :s perpendicular to the Sfow direction and the spacino in the ?est direction is d. Pref- 
erable spacing Is selected in the slow direction In such a way that interpolating between two adjacent sample points A 
and B m a singie line parallel lo %\® slow direction would have an enor rx? greater than the error Incurred by interpolate 
ing between adjacent sample points In any of the merged ordinal functions. The specif 5g In the slow direction should 
therefere be such mat If a line parage! to the fast direction of any of the merged original f unclbos Is used as a coordinate 
ax?s, then the distance along tiis ooordin:ate ads between A and B wiff be no more than d. If san he seen from Fig. 4 
that ttxr desired distance between ad^acem polnls in the slow direction of f | is therefore 

ja\ 

Fig, 4 illuetrates choice of spacing m the slow direction The angular diflerence in radians between the slow direc- 
tion (verticai lines} of the function f |. , ounentty bf^. formal, and the slow direction (solid diagonal lines) of any of the 
merged ojlginal function Is never neater t?ar^ ! ~i)xf{2Q} The fast direction of the j^erged original funotlc^ Is pen 
pendlcu^ar to Its slow one and its .spacing- in that direo^on Is d. The distance In that direction between adjacer?t san^p^e 
points A and 8 of / £ sbonld not e^eed d. On the left, the -first merge l>?vel is shown On the right, an arbitrary n^erge 
levei Its shown. 

Since the diameter of the support of the- reconstruction is Nd f the number .of sample points necessary along « 
gle line in the slow direction of ? £ Is Nd divided by the sampling if nerval or 

N s:n( (2 ; - ? ). 
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Let t*s that if Ihfe opo-sber Is not an integer then the lowest integer greater than tNs nun^er will bs used. 

Since the sample spacing, in the last direction & a&*$y% d, \tk$km& &m toi any I and J, there are H lines parallel to the 
slow direction alon« whioh f i w$ be sampled, so she total number of points in gj & 



N f N sin (( ?J " l)rrf ( 20) I < 



- • *" ■ 

Ate, there are G/2 : different grids 9 J- to be co^p^ecs (i.e.. a! the ith level of merges, 3 takes ih« values 1 . .,G/£ : ). It 
Worn* ft<$n Ihis that the total number of poifll. values io h®-GOWpi$ed- and stored at the Ith level of marges is 



The sum of fnls for all levsis of merges, La. for i » l^.Jog^G, us less than 

» /V i! tog 2 Q + /vD> 

If, as is normally assumed, O * R then the order ollhe bac4<projeeti6rj is hj-lbgN. 

£i> During Implementation of t&s sl§or?lhrn some pr&c&sal Issues should he considered. For i =* l % . \®Q>Q, samples 
of f f are ooarsput^i^W: g/); ? and g g } by ihler poi&lion< but there will ^ehfcrally be some points where /{ i&to be sam- 
pled that era not sunoorsded t>y sample points of f ,;Vl or of f %) . and extr&|x>lation mus* be performed -there Instead of 
Interpotei&n. Computing a few sanies of f\ just outside tie di£K In which the image is to be reconstructed will reduce 
the nunrtosr of points at which extrapolation & accessary at the {\ * IHh level, but wHI In itself require eatrapoletioo 

so unless a sufficient .number of samples at points ou&id® the disk were computed at the (M) -th leva? A with the 

same computational complexity Is conceivable in which at each iev«l of grid merging , sample points outside the disk' are 
chosen in such a way that exmoolstioo Is never required 
in me present tniplsr>w&atk5o.. one mors than the 



points per line m the slow direction p:eeorit>ed above is ased. ao tiat the first and iasf sa^i?^ posnfe an asoh §ne 
would he at or beyond the boundary of the dek supporting the ra<:onsfn.sct?on. Ttiia does not change tf^e cnn^^af ionaJ 
co?npiaxity. 

W« have found empirically that jrnage quality can be improved by dousing the ntinte of samples pe* Radon pro- 
;eet!on by Intarpol^lon, ssr by doiMn§ the nunt^er of pr<sie<^k>ns (i.e., the mimbdr of angles at tth£h pmi^iona are 
coripoted) by Interpolatlor?. A more ef ioieht way to in-prove image quality is to sempte f{ lor alt I and j st no fewer than 

*$ S points in the slow Predion, although the expression for the desired distance between adieceni poir^ts would Irr^y that 
fewer ean-ple points are necessary in the slow direction in the fi?^t few levels of rnerpee. In any case., neither of these 
ao^ustmems ere neooss^ry when the Postprocessing correction n-ethod described l^eio^ Is applied. 

Images rosuill?^ the mu&scale beckproiection slgoh^?m as described up to this point ere somewhas blurred 
due to the many lrr*erpola^ons neoossary The point spread function oi the slgonthm wider than that of the classical 

§q bac$^ro|eotso?t. We have found that an G{N : '*) Fourier domain correction of the Image greasy reduces the width of the 
poir^t spread function of Ifie multiscale beokpi^ectlon and enhances the resulting images A similar correction can be 
performed after the classical hackpro^eclsom but with less of an imprpyeniem> and images procsuced by the rnuitlscale 
method with Fourier don5a?ri correction are aileas<t as good qualitatively as those produced u&lng rhe classical bac^fo- 
j action with Fourier .domain coo*eotion. 

m If the po?rn ^ e$sd 1 unate or a given dacHprojeclion method were shift Invariant, then the obvious correction would 
be to divide the Fourier ^ansform of me reoons^L^ction by that of the point spread lunction and to take the Inverse Fou- 
rier transform of the result However, the point spread Junctions of both the classical end muitlscale backprolectlon 
methods vary sightly over the inuxgo. Postproc^sing correc^on herein coasts of dlvxilng the F$>urler transform of the 
reconstruction by a Gaussian whloh appmxln*&ies the Fourier iransfcrms of the point spread furscttdns obtasned al var- 
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sous points In the image, The division k> p^onned in a disk around the origin of tha Fourier plane. Outside this disk. 

for higher frequencies., the Fourier transform of the reconstruction :a retaioad* and no dfrfeton Is performed The 
width of the* Gaussian & chosen in suck a way as to optlm&e tftfe approxm>sik>a As we mil show, this technique Is more 
elective for &&myHtea&i$ bae^xojectlon than for the e&ssicsi one as the po^l: spread functions obtained for the muh 
tlscak* method can be mors dotaly approximated by a Gaulish. 

For each of tha two methods under rans*dera#on {the ctet$$&a! and muStisoaie backprojeotion methods}., the 2 D 
Gaussian was found wh&h fit a sanction of point responses better than any atft*r£~0 Gasman. Mora precisely the 2- 
D Gauss^ans are o? me form 



where the wfcfth o>0 of theGaussjah may bo chosen arbitrarily. In the context of the present Mention, x and y may be 
regarded as »M c&orc£ngi&$ the x coo:€fcoaies of horteon&gy adjacent pteefs differ by % and the y coordinates of 
vertically adjacent plxe& differ by 1 },. and then the goal was to choose trie vaJuo o 0 of o for whioh the resulting Gaussian 
best fit a selection of point responses in the following -sense 

The point responses wars all rapf$sentaci by re&f -valued matrices. These matrices were averaged, and the central 
? x 7 submatrix A,, with the peaks of the point responses at its center; was extracted and normailsad. Let the Indices \ 
and j of A run h-om -3 to 3. Let 



* </ r 



For each method, o^ Is dehe^d as the vsalue ot o which minsm&es 



Let us assume that baoKpmjec&on is performed on a square grid of pixeis, or on a circular oootr&l region of this 
square The computation of the Radon projections of trie image contained in this grki and tha convolution arxs backpn> 
joeibn oi the Radon projections bask on to the grid, can be performed in sud^ a way that the entire sequence of oper- 
ations is symmetric with regard to the hcvlsootaf and vortical axes of tha orid and its two diagonals, in such a case., a 
point response anywhere m the grid is Identical, up to rotation end rehee&an, to as many as coven other point 
responses, whose positions are chained from the position otlhe original point response by rejection across one or 
more ot the axes of symmetry of tha operates. See Fig. 5. As shown in Fig. 5, a point response of tha Radon trans- 
fbrn>CB sequence is Identical, up to rotation and ref action, to as many as seven other point response It follows that 
whan selecting positions on me grcl at which to compute point responses kit me purpose of computing <s r> one may 
restrict the positior^s to a single octant of the grid and compute the point responses si those positions, and than tie 
point responses st the corresponding positions in the other seven octants can be obtained by. rotating and reflecting 
confuted point responses. 

In our research, point responses were computed at 1-5 randomly &el&oted points in a circle octant and the point 
responses at the scw&eponding points in the ramainsng octants were obtained by rotation and ratiactte of the original 
I S point responses, as oxpiained apc^e. 

Point responses and the r«suH£ng value of o^ da^end on the process by whicn me Bedon prpjoctiona are stained, 
in oor research, Radon pn^eotlon datB were computed tr\ software, but an apparatus in^^meming tha m^tpd 
?:less?:ibad here mioht obtain the projector; data from physical measurements. In any case for any suct^ apparatus. o, v 
should be cos-nputed from point responses generated by Radon data obtained from the source from which the appara- 
tus will obtain the Rsdon data In practice. 

Fig, -6 indlca^ee the optimal vaiue of <y t arxj the decree to ^tiich the resulting Gaussian fits t^>e point responsee tor 
the csssslcal backproj actions and several variants or hie mu^soa^ backprojsctlon. Over each point response grayscale 
Imago Is displayed a pa^r ot numbers. The first Is the optima; ^alue o c of o lor tha method, and- the secor^ Is the error 
i: ^Iv^n by 
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II ca^s. be r*>an ffe&f w» o% ; is gre&ter tor the muK&c&te bat^xojactlon than for the dasslcaj baokprojeetlon, Implying 
that rnt£tte$Gtl& baewro*adfi«o without posiproces$lr*g has a wirier point response and "thersfom produces tmxner 
imaga. the rrtuft&cs&a method point rr^per^as can be approxi m&tad betf&r by a Gaussian, and therefore she nu&lscaie 
bad^5ro|e^fe *>lih ps^l^oo^^ng can produos sharps* images than m& clashes! baefcprojectos with po^M&s&ng. 

Rg. ? stas sun&oe pio& of point spread functions fear the tia&s&ai End muHiscaie bac^ro^ofcbn methods and 
surfeo* plots or ifc**r Fourier transforms, ft can be seen from these figures- that although iba uixw^ed rm&iscaie 
bac^roiaotirx; product a wfcter po*nt apread than the <tesafc&I method, the Fourier transform of ttte rm^scai* point: 
spread function mora cioeety rsssmbfes a Gaussian than the Fourier grdn&ferm of the classical point spread Unction. H 
foitews that the mu&sseaie mathod better fends itself to the postprocessing correction described bare than does the das- 
stc&i method. 

The corrddksn lust deser&sd as a postprocessing con-action can stoso be erfaeted during the convoiufion step 
instead of requmng a post-processing st«p> based on the toflowirQ. r|- Allows easily from the Radon Stee Thaorem (2) 



that 



f*f{p) :: >/T^f( p cos g-> pjm 0) 



Thus the pootbeK^prcjao^ procedure of dividing by the ..&»q discrete HXttksr transform of 

'> 

in tha disk in ?be Fourier domain defined by 4- y e ' < fr may be replaced by division at 
by the 1 -D discrete Fourier transform of 

in thsirhrery*! defined by p 2 < R* This method of applying the correction may be preferred "m many cases -as ft does 
not require introduction of an eotira add&ona* step and also does not require the amounts of data storage oacassary 
for a 2-D discrete Fbuher transform. 

It is import&ntto note that the computations of <y Ct> of tha S-D discrete Pawnor transform of 

> 

and at tha t-D discrete Fourier transform of 



EF0 811 052 A2 



can be dons ,r ofHlne K r that ie ; they can bi? done before sensor data are read in for the purpose of nxconsfn^otioo, aid 
therefore it is only the division of the 2-D discrete Fourier transform of the Imsage by the 2 D discrete Foun«f transform o! 



or the dlvlsson of 

p#f(p) 



by the 1-D discs^le rouner transform of 



that by necessity add to the amount of time tie method will require to effect the reconstruction. 

Ha 8 shows an image from uvHcfc Radon data ware computed -for the purpose of testing efferent variation of OB. 
Fig. 9 shows the resets of OB with classical backproject&ns. Fig, 10 shows the resume of OB with mufecale backpay- 
lection without postprocessing &itwirft doubling of the date by means of Interpolation. Fig . 1 2 shows the results 1 of CB 
with mulflscaie b&ckprp^ectlon, without doubling of the data hut wrf vpostpK^essing. 

Although the invention has been slxm j n and described *n terms of a specific preferred embodiment and variants, 
charges arsd rrxcfcflcatlons are possible which do not depart from the spirjt, scope, or sontemp&tion of the inventive 
concepts disclosed and taught heresiL Such are deemed to tali wMm the purview of the inversion as claimed. 

Claims 

f v. A method of computerized torr^raphy, comprising the steps of: 

(a) scanning a subject by projecting radiation toward the subject ; 

(b) sensing the. projected radiation with, a plurality of sensors: 

(c) prooesr^ng the output ol the sensors to oreale a matrix of samples or an Image by: 
«!}• § ilterin§ the sanies of each Radon projection sample vector ry where I « 1 ^G, by: 

i. computing the discrete Fourier transform r of r 5 using an algorithm of ord^r N log hi where N is the length 
of the vector: 

11 multiplying the jfh element of r * by j, for j < N/2> and multiplying the Jth element of r 5 by M-j, fcr j > W2, 
where the elements of r* are numbered 0 through M-i : and 

lit. cottpuitag the inverse discrete Pouher transform g c J of the modified r* using an algorithm of order N log 
and 

(e) performing the k1h level of merges, for h « 1.,., ,.log>A by compiling the grid gf t for i « 1 by m&g&s 

arKt by means of rotation of coordinates, interpolation, and addition. 

2. The method of claim 1 therein the algorithm to compute the discrete Fouber transform Is the FTT algorithm and 
the algorithm to compute the Inverse discrete Fourier transfor m Is the I.FPT algonthm. 

3. T he method of claim 1 wherein the number of samples In each Badon projection vector r* is increased by means of 
Interpolation, -and/or the number of Radon projection vectors r 5 la increased by means of imerpolallon, before fie 
st$p of computing the discrete Power transform Is computed, 

4. The method of claim t % wherein the number of samples, of each vector g^ Is increased by means of interpolation, 
and/or the number of vectors §f Is Increased by means of inter poiation, after the step of confuting the Inverse dis- 
crete Fourier transform Is performed. 
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& The metnod of ofasm 1 wherein the veclo? r j Is mii^l&d by a weighting function «o the st>sp 0? {e} n, 
a, A method of computerised •omogmphy comprising Ihe steps of; 

5 (a) seaming a subject by projecting radiation toward the subject : 

(b) . sensing the projected radiation with a plurality of sensors; 

(c) processing the outpiri data of the &&r&ctr& ?o sreafe a matrix of samples of m ins&ge by: 

(d) convolution and hacl<proj<&ci*an of the cnjttput: data; 

re) computing a selection of the= point spread functions produced by the Radon tfjxnsfoj-?rvoonvolutJO?vbac^o* 
v>> section suite; 

(i) reputing the width ^ of the Gaussian that most oSoeefy tits the computed pc»m spread functions according 
to « matching mathematical edtehon; 

(0) ^npufing G , me 2-0 derate Fourier transform of fn» M x N of samples of the Gaussian of width o 0 using 
an algorithm of order no greater than N^iogN; 
*s (h) computing the 2-D discrete Fouher transform of the N x N image matrix using an algorithm of order hriogN; 

(h d»v$c$n§ each component of the resulting matnx by the corresponding sample of G> with treatment ol corn- 
portals at which G is near saro; 

(j ) compiling the inverse 2-0 discrete Fourier transform of the resulting matrix, 
(k) ebspi>ayln§ the results on a visual display 

7, The n^ethod of cteam 8 wherein steps, (h) through (j) are omitted arid the function 
01 me convolution & dMded by tie discrete Fourier transform of 

x y cr , 

lor all p such that p 2 s for soma R 

35 0. A. computerised tomography scanner, comprising: 

means tor soenning a subject by projection thereof wsth a source of radiation : 

$b) means for sensing the projected radiation wish a plurality of sensors on the side of she subject opposite the 
40 source -erf radiation: 

<c) means lor processing the output of the sensors to create a matrix of samples of an image, including 

(rf> means for computing the discrete Fourier transform ? ol t l using any algorithm of order N tog N. where Mis 
4$ ihe iength of the vector, 

ie) mesm 'for multiplying jth element of r ! by j. for j < and rriultiplyir^ tha Jth element of r by M ■ i> for j > 
M& % where the elements of r are numbered 0 through N - 1; 

m (1) means for computing the inverse discrete Fourier transform of the modified ? using arty s&otifcrc of 

o:'der M sog 

($} means lof cor^Uting ^ g^-fer Z * t ,Q/2 k . by merging gg:\ m&§%j by means of rotation of coordi- 

nates> Interpcsatlon, .a?x! add-on; and 

<:c: 

#\) means for creating a viavsai display of the matrix of sampies. 

T ne a^^rete of cialm S ? wherein the means So? conpu^r^g the discrete Fourier transform is me FFT aigorithm and 
the means for computing the inverse discrete Fourier transform is the iFFT algorithm. 
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10, Tii* apparatus of claim 8, wherein the number of samples In each Radon projection vector r Is ^ms^ by means 
of snt^poLatioo ; aod/dr the number -of Radon ptagdcSon vectors r is increased by m«§ of Interpolation, before 
computing fhe discrete Fourier transform. 

s il. The appa? atus ot o\mm B r wherein the number of samples of each vector g*j is increased by means o? fote^xjS&tsosx 
and/or the number of vectors l§ irK;re&sad by means of inte? po^en > after computing the inverse discrete- Fou- 
rier transform. 

12, Tha ^.paratue of $, wt»*c«ta the vector r is multiplied by a weighting function in m<sane to rmj&iplyihg. 

m 

13. A conpuierl-rad tomography scanner comprising: 

(a) means for scanning a subject by projection &t radiation: 
ia (b) means for sensing the projected radiation mih a plurality of sensors;. 

(c) means lor processing the output data of the sensors to create a matrix, of sanies of an image , Including; 

(d) means *or computing the convdutlon and backprojection of the output data; 

(a) means lor computing a selection of the point spread functus produced by the Radon tmnsform^oonvolu- 
fiombackpmjecfion suits; 

(!) means tor composing the wsdth of the Gaussian that most cio&efy fits the combed point spread functions 
2$ according to a mathematical criterion: 

(g) means for confuting Ci, the 2-0 discrete Fourier transform f2DDF*0 of the N x M matrix of samples of the 
Gaussian of width a.,; 

m (h) means for computing the ?>D discrate Fourier transform of tie N xN image masri* using any algorithm of 

order hhlogN; 

{1} means for dividing each component of the resulting math* by the corresponding sample or G, with treatment 
of consonants at which 3 is near sere: 

rj) means for computing the inverse 2 P discrete Fourier transform of the resulting matrix: and 

fk) means tor converting the matrix: data to a visual display 

m 14, A computer readable memory medium encoded v*&* data re^rese^ing a ceroputor orograro to use wrm a cornpm 
termed tocography scanner and a computer 1o gan^ate a visual image by:: 

ra) fiitenng samples of Radon projection sampl e vectors r\ where I » !,,...Q. by: 

46" ?: computing me discrate Fourier transform r' of r 5 using an algorithm of order N log N> where hi Is the lengfh 

ol the vector; 

H. motiving the |th element of r by j, for J < N/2, and mu&plymo the jth element ot r ; by N-J, tor i > 
where the element of r are numbered 0 through N-i ; and 

iii. computing the inverse discrete Fourier transform g9 ot the mcMm r : using any algorithm of order H 
so log N: and 

(l>) pehorming the kth level of mergas, for k «< 1,...jog^G\ by computing the grid. gf . for I t v . : .»<S&*> by merging 
q£;\ and o^l hv means of rotation of coordimtes, interpolation, and addition, 

5^ t& Toe computer ratable meinory medium ot claim 14 sn the- form of a floppy disk, an ASIC, compute? chip or circuit. 

16, The computer readable memory medium according to claim 14 wherein the number of samples In eaci^ Radon pro- 
taction venter r 5 :s increased by means of interpolation, and/or the number of Badon proieclkin vectors r s ' Is 
iOCfeases^ by means ot Interpoiafion, before tt^@-^ep of computing discrete Fourier trsralsxrm is c<snpu£ed : 
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IT. The con^xster 8$&d§bfe memory medium according jo cteant 14 wherein the number or samples of each vector 
is increased by rmans of bierpolatiom and/or the 'number of vectors g^j io incased by means of i*terpo^ors ? 
after the step of oonptissrfg invert discrete Fourier transform Is performed. 

1& The computer reatkbte memory medium.. according to claim 14 wh«rdn Ihe vector r- :s mulled by a weighting 
ksnctioo in the step of rcu&p$yin§ discrete Fourier transform is performed. 

The compter resdshte memory medium encoded vvith data representing a odf*puter program thai can cause a 
conipuler to function so execute the method of: 

(a) computing a sek-o&on of she point spread functions produced by the Radon h^-^form-convoiut5on-backpro- 
jectioo suste; 

(b) computing the wklth o^ of th& Gaussi&n that best nts the confuted point spread functions according to a 
m&rrtem&fccas cntenon; 

(c) computing G, the 2-D discrete Fourier transform (2DDFT) of the H x N matrix of samples of the Gaus&san 
of width <y c ; 

(d) computing the SO&FT of the N x N image matrix using an algorithm of order N 2 iogN; 

(e) dividing each component of the resulting matrix by the corresponding sample of Q c wth apeciai treatment 
of components at which G is nee? zero; and 

(!) computing the inverse 2 D discrete Fourier transform of the resuming matrix. 
20, Tho eoe^puter readable memory medium of claim 14 in the form of a floppy disk, -an .ASK* confute? chip or orculi 
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